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TECHNICAL NOTE 4277 

A BODY MODIFICATION TO REDUCE DRAG DUE TO WEDGE 
ANGLE OF WING WITH UNSWEPT TRAILING EDGE 
By Wi lli am C. Pitts and Jack N. Nielsen 

SUMMARY 


Ward* s s lender-body- theory formula for zero-lift drag contains three 
integrals plus a base-drag term. Two of these integral terms depend only 
upon the cross-sectional area distribution of the body. The third inte- 
gral term depends only upon the body shape and axial slopes at the base 
of the body. This term is neglected in the transonic area rule because 
in many cases it is zero; however, there are also many cases in which it 
is not zero. This paper examines the term for the possibility of drag 
reduction for a particular case. The model considered consists of a 
body of revolution in combination with any wing that has an unswept trail- 
ing edge and a constant trailing- edge angle along its span. It is found 
that (neglecting any change in base drag) a drag reduction is obtainable 
which, for the case considered, is an additional 12 percent of that 
obtained with the area- rule modification. The probable effect of viscosity 
on this theoretical result is discussed. 


INTRODUCTION 


The transonic area rule (ref. 1) relates the drag of a configuration 
qn rt the drag of an equivalent body of revolution having the same area 
distribution. Engineering methods have been developed from the area rule 
for calculating the drag of airplanes and missiles at zero angle of attack. 
The basis of the area rule in slender-body theory can be investigated by 
studying Ward’s drag formula (ref. 2). It is found that the equivalent- 
body concept holds rigorously only if certain conditions at the base of 
the body are met, and if the trailing edge of the wing is swept or cusped. 
Frequently these conditions are violated, as pointed out in references 3 
andi 4, and additional drag is obtained above that of the equivalent body. 
Beradt (ref. 3) states that two bodies have the same drag only if in addi- 
tion to being equivalent in the sense of the area rule they have the same 
cross-sectional contour at the base and the same streamwise slope around 
that contour. In reference 3, Bemdt investigates how the difference in 
drag between equivalent bodies depends upon the difference in base shape. 
T>»p important conclusion is that, within the limitations of his approxi- 
mate theory, the drag of a slender body having a finite cross-sectional 
slope at the base may he considerably reduced by spreading out the base 



2 


NACA TN 4277 


contour from a circle without changing the distribution of cross-sectional 
area. Lighthill (ref. 4) has evaluated the drag increment over that of 
the equivalent body of revolution for planar and cruciform wings with 
unswept trailing edges alone and in combination with cylindrical body. 

It is clear from references 3 anu 4 that the body shape given by the 
transonic area rule does not give the minimum possible theoretical drag 
for all configurations. An example is a wing-body combination for which 
the wing trailing edge is uncusped and lies in the plane of the body base. 
The purpose of the present paper is to determine (within the accuracy of 
Ward's drag formula) how much the drag of this wing- body combination can 
be reduced by modifying the streamwise slopes, of the body at the base. 
Since viscosity can have an important effect on the reality of the 
inviscid-fluid-theory results presented herein, the probable effect of 
viscosity is discussed. 


SYMBOLS 


a 

body radius at 

X = l 

a o(x) 
t>o(x) J 

parameters In 


a 2 (x), \ 

a 4 (x), ...J 

parameters in 

cpja expansion 


b 

c 

C P 

D 

AD 


wing semi span at x = l 
wing root chord 
pressure coefficient 
drag 

drag after adding shape modification minus drag before 
adding shape modification 


*o(x), 
d 2 (x) , 

fo( x )> 
fa(x) , 

ga(x), 

S4(x) } 



parameters in cp^ expansion 


parameters in cpj_ expansion 


slope- amplitude function for shape modification 


l 


length of wing- body combination 
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m,n 

q. 

M 

R(x) 

R(x,0) 

s 

S(x) 

V 

X,r,e 

Xl 

X 

Z(x,y) 

I 

9 

Sep 

Sn 

A 

B 

C 

i 

m 

W 


positive integers 
free- stream dynamic pressure 
real part of complex function 
surface of Basic "body 

surface of Body with, shape modification 
length, of arc along contour 
cross-sectional area of model at x 
free-stream velocity 

cylindrical coordinates (see sketch (B)) 

axial distance from leading edge of wing root 

complex number in x = l plane 

upper wing surface 

dummy variable of integration 

perturbation velocity potential 

outward velocity component normal to surface 

Subscripts 

potential after adding shape modification 
basic body plus area-rule modification 
basic-model combination 

interference between wing and cylinder of radius a 

shape modification 

wing 
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Superscripts 


first derivative with, respect to x 
second derivative with respect to x 


PHYSICAL CONSIDERATIONS 


The general statements of the Introduction can he given specific 
meaning if the drag formula of Ward is considered. 


D 

q. 



S"(x)S ,, (i)d| dx + 


- Cp s( z) 

Base 




The first two terms depend only on the distribution of cross-sectional 
area along the length of the wing-body combination. The third term 
depends only on the. shape of the wing-body combination in the crossflow 
plane of the trailing edge (contour C in sketch (a)) and the streamwise 



slopes of the wing and body surfaces approaching the trailing edge. The 
last term is the base drag which is not considered, although some of the 
body shape changes considered might possibly induce significant base- 
pressure changes. 
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Since the first two terms of equation (l) depend only on the area 
distribution, they represent the drag of the equivalent body. Frequently, 
however, only the first term is used in engineering methods for calculat- 
ing wave drag based on the equivalent body concept. It is, therefore, 
important to know when all terms (other than base drag) but the first are 
zero. Two cases can readily be found. The first case is that of a body 
with a pointed base, or a wing with a cusped or -swept trail ing edge. 

Then S’ ( Z) is zero and the contour C is zero. The second case occurs 
when the body and/or wing is tangent to the cylindrical extension of the 
contour C. In this case the second term is zero because S’ ( Z) =0 «.t>h 
the third term is zero because Scp/Sn = 0. For other configurations, 
the second and third terms usually contribute to the drag. 

The purpose of this paper is to investigate what drag savings are 
theoretically possible through control of the third term. This is done 
by modifying the streamwise slopes of the body on contour C wi tho ut 
changing the body cross-sectional area. The streamwise slopes forward 
of x = Z can be chosen arbitrarily to fair into those at x = Z with- 
out affecting the drag. The first two terms of equation (l) are unaffected. 


ANALYSIS 

Velocity Potentials 


The basic model and coordinate system used in the analysis are shown 
in sketch (b) . Only the shape and slopes of the model in the x = Z 



plane need be specified. The model shape forward of x = Z is arbitrary 
except that the model must be slender in the sense of Ward’s theory. r T t hp 
area-rule modification is included in the basic model. To simplify the 
analysis, the case in which the wing- trailing- edge angle is a finite 
constant along the span is considered. 
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The perturbation velocity potential in the crossflow plane of the 
basic model can be written in the form 

<P C * ‘Pw + ^ + 9 i ( 2 ) 

where qfy is the potential of the wing alone (including the portion 
blanketed by the body) , cpg is the potential of the body plus the area- 
rule modification, and cpj_ is the interference potential that cancels 
wing-alone components of velocity through the body surface. To cp^ there 
is added a shape modification potential, cp^. As previously discussed, 
the purpose for the shape modification is to modify the streamwise slopes 
of the body at x = l in such a manner as to reduce the drag of the 
wing-body combination. For a reason which is subsequently pointed out, 
the restriction is placed on the shape modification that the body cross 
section at x = l must be a circle. This is not a serious restriction 
because any meridian slope distribution can be faired into a circular 
base by properly shaping the body forward of the base. The potential, cf^, 
can be added directly to cpp provided it does not violate the boundary 
condition of no flow through the wing and body surfaces. 

The wing-alone velocity potential in the crossflow plane, x a l, is 


m . • I Z« ( Z)m[(X+b)log(X+b)-(X-b)log(X-b)-2b] (3) 

W 3t 

where X = re^® is the complex variable in the crossflow plane. The 
other potential components, q}g, cp^, and q^, are obtained from the general 
slender-body potential 


<P 


00 

ao(x)log r + b Q (x) 

m=i 


cos m0 


( 4 ) 


This series converges for r greater than the body radius. This is the 
exact slender-body theory potential only for body cross sections which 
are circular. For this reason the restriction is put on the body shape 
modification that the base of the body shall remain circular. 

The parameters ao(x) and b Q (x) are functions only of the model 
cross-sectional area distribution. The parameters a m (x) depend only 
upon the streamwise slope distribution of the body. Therefore, the 
potential of the body plus the area- rule modification is of the form 


q> B = a 0 (x)log r + b 0 (x) 


( 5 ) 


and the shape modification potential is of the form 
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<Rm 


00 

-I 

m=i 


^mOc) 


cos m B 


( 6 ) 


The parameters a^x) are determined from the boundary condition 


Then, 


00 

T 

m=i 


VSr/ r=a 

mamCx) 


= VR*(x,0) 


□hh-x 


cos m0 = VB’ (x,0) 


This suggests a shape modification of the form 


R» (x,0) 


00 

-I 


Sm,(x) cos m0 


m=i 


where g^(x) is a slope amplitude function. To the order of accuracy of 
Ward’s formula ( eq_. (l)) this type modification meets the requirement of 
no body cross-sectio n a l area change, and if only even values of m are 
taken, it satisfies the boundary condition of no flow through the wing 
surface. Even values of m are also required for symmetrical flow about 
the (x-z) plane. Then, 

V a saC x ) _ „ . , 

^ = )_ J ~^— COs2xie (7) 

n=i 


with 


a 2n( x ) = -V a 2n g2n( x ) 

The normal velocities produced on contour C by the n = 1 arnS n = 2 
components of cf^ are shown qualitatively by the arrows in sketch (c) . 



Sketch (c) 
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The iirberferen.ce potential cp^ is obtained from cp^ "by using the 
boundary condition that there be no net flow through the body surface. 
Mathematically this is written 



( 8 ) 


Expansion of equation (3) in a Fourier cosine series that converges for 
r < b and 0 < 0 < if gives 




CO 

= d Q ( l) + ^ danC 2) cos 2n0 


n=i 


vhere 


do ( 2 ) = ^iog b + | - 1 ) z* ( 1) 


cLsnO) 


r>2n 


4 r 


Z*(2) 


(9) 


* Ln(2n-1) b 2 * 1 4 n 2 -l 1 . 

Then cp^ is obtained from equations (4), (8), and. ( 9 ) in the form 


9i 


OO 

* foO) lo S r + Yj 


fan( 1 ) 


-an 


cos 2n0 


( 10 ) 


n=sx 


where 


foO) 


— Z* ( 2) 

Tt 


f2n( l) 


a anfi v / aga^x _2_\ , , » 

(an-l^L 2 *" 1 ’ 2n+V 


Theoretical Drag Reduction Due to Body Shape Modification 


As previously discussed, only the third term in equation (l) enters 
in the drag increment produced by adding a body shape modification with- 
out altering the cross-sectional area distribution. If subscript C 
refers to the potential of the combination before adding body modifica- 
tions and subscript A refers to the potential after adding the modifi- 
cations 
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f< p, ^ as + ^ f 

4 v 2 Jn A Bn V 2 X 




Sqpg 

Bn 


ds 


(H) 


■where 


Cp A = CD- + Cp 

\A m 3 Xm 


(The contour C is sho-wn in sketch (a).) With the conditions that 
_ a<PB acpi ^PW Bcpi 


Bn 

body 

Bn 

Sn 

AD 

2 

r* 


Sn 


^ + ( V'Pb+'Pi+'H.) % 


a de 


Bn ~ Bn 
4Z» ( z) ^‘ b 


on the 


/ V* 


( 12 ) 


Insertion of equations ( 5 ) ^ (7) ) (9); and (10) and the orthogonality 
property of the cosine function into equation (12) gives 


q. 


= P Z X [ dan ^^ f2n ^ a2n(z) 

n=x 0 


cos a 2n0 d9 - 

r 2n+i 


4z* ( z) y r 10 Wz) 

V LJ a r 2 * 
n=x 


Integration gives 


AD 

q. 


a 2 


00 



n=x 


«sn(2)Z , (l) 



(13) 


We obtain the optimum values of &sn. ( l) hy considering each term of the 
series separately. The values of gsn(l) that give the maximum AD/q. 
are 


82 n( 7>) = “ 


4 

(2n-l) jt 



&2n -A 

h^-V 


Z’(Z) 


(3*) 


The maximum drag reduction is then 
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AD 

qa 2 


2 

ir 



n=i 


2 

2n-l 



a 2 n- IV 
^an-iy 


|2 


( 15 ) 


As expected, the optimum value of g 2 n( Z) and the drag reduction are zero 
when the wing trailing edge is cusped. 

Although the total drag increment depends only on conditions at the 
model "base, the drag increment does not originate there. Actually the 
drag increment is distributed over the entire winged portion of the model 
as illustrated by the following example; Since the distribution of the 
amplitude function, g 2 n(x), is arbitrary for x less than Z, choose 

ge(xi) = Ai + A 2 (%y 


where xi ;= x-(2-c) a; ad Ai and are chosen so that g^( 2) satisfies 
the condition of maxi mum drag reduction (eq. (l4J). The distribution of 
drag increment due to the addition of this g 2 ( x ) shape modification is 
shown in sketch (d) for a biconvex- section rectangular wing. 




Sketch (d) 
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Inf luence of Boundary-Layer Shock-Wave Interaction 
on Drag Reduction 



Thp foregoing calculations are "based on inviscid-flow theory which 
ignores the interaction between the trailing-edge shock wave of the wing 
(and "body) and the boundary layer approaching the base. The questions 
arise as to how this interaction affects 
the over- all drag and how it affects 
the drag reductions due to the type of 
modifications of the body shapes cal- 
culated herein. In the theory the 
streamlines at the trailing edge are 
assumed to be as shown at the top of 
sketch (e). The streamlines are par- 
allel to the sides of the trailing- 
edge wedge and undergo a pressure rise 
in traversing the trailing-edge shock 
wave. Because the pressure rise 
occurs behind the wing, it does not 
act to decrease the drag. If the 
trailing edge were cusped, as shown at 
the middle of sketch (e) , the pressure 
rise would occur through a gradual, 
compression in front of the wing 
trailing edge and the drag due to the 
trailing-edge angle would be 
eliminated. 


Inviscid S' (2)*0 



Inviscid S' (2) s 0 


The influence of viscosity on the 
drag is somewhat similar to that of 
cusping the trailing edge. The bound- 
ary layer allows the pressure rise 
through the trailing shock waves to be 
transmitted upstream. This thickens 
the boundary layer; compression begins 
over the wing surface; and shocks move 
up in front of the trailing edge. The 
accompanying pressure increase over 
the rear of the wing acts to decrease 
the pressure drag below its value on 
the basis of inviscid fluid theory. Such an effect has been shown by 
pressure distributions (e.g, , ref. 5) for two-dimensional airfoils in 
supersonic flow. The amount the pressure drag is reduced below the 
inviscid value depends on the boundary- layer thickness approaching the 
trailing edge, the shock- wave strength (Mach number) , and whether the 
boundary layer is laminar or turbulent . A phenomenon similar to that 
shown at the bottom of sketch (e) is also important in the wave-making 
resistance of boat hulls. As pointed out by Havelock (ref. 6), the 
wave-making resistance is reduced by virtue of a thickening of the 
n frietion belt” near the stem. 



Sketch (e) 
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Since "boundary- layer shock-wave interaction influences the flow at 
the trailing edge and since the drag increment expression (eq. (15) ) 
depends only on conditions at the trailing edge, it is pertinent to ask 
if the drag reductions computed herein on the "basis of inviscid fluid 
theory are realistic. Although AD/q is a function only of trailing- 
edge conditions , the drag reduction is distributed over the entire winged 
portion of the configuration as shown by the example in sketch (d) . For 
this example the total drag increment up to x x /c = 0 .65 (shaded region) 
is zero, so that the net drag reduction given by the present inviscid 
fluid theory is distributed in the region 0 .65 < (x x /c) < 1 . Since 
boundary- layer shock-wave interaction effects will influence this drag 
distribution only near the trailing edge, it appears that most of the 
drag reduction predicted by the inviscid theory can be realized; the 
actual amount can only be determined by experiment. 


CONCLUDING REMARKS 


Equation (15) gives the maximum drag reduction obtained (in inviscid 
theory) by modifying the streamwise slopes of the body without changing 
the cross-sectional area. This expression is independent of the wing 
plan form since it depends only upon conditions at the base. However, 
the drag reduction is distributed over the winged portion of the wing- 
body combination to the extent that the streamwise slopes of the body are 
modified. This drag distribution does depend upon wing plan form. 

To give an idea of the order of magnitude of the drag reduction 
given by equation (15) , a comparison is made with the drag reduction 
given by the area-mile modification. For a delta-wing cylindrical-body 
combination (sketch (a)) with b/a = 3* "the additional drag reduction 
(neglecting any change in base drag) given by the g 2 ( l) shape modifica- 
tion is 12 percent of the drag given by the first term in equation (1) 
(the area-rule term) . This percentage varies somewhat with different 
configurations, but it should remain the same order of ma gni tude. The 
effect of the other Ssn( 1 ) modifications is negligible compared to 
the g 2 ( l) modification. 


Ames Aeronautical Laboratory 

National Advisory Committee for Aeronautics 
Moffett Field, Calif., Mar. 21 , 1958 
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